Frequency doubling and memory effects in the Spin Hall Effect 
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We predict that when an alternating voltage is applied to a semiconducting system with inho- 
mogeneous electron density in the direction perpendicular to main current flow, the spin Hall effect 
results in a transverse voltage containing a double-frequency component. We also demonstrate that 
there is a phase shift between applied and transverse voltage oscillations, related to the general 
memristive behavior of semiconductor spintronic systems. A different method to achieve frequency 
doubling based on the inverse spin Hall effect is also discussed. 
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In optics, frequency doubling (also called second har- 
monic generation) is obtained from nonlinear processes, 
in which the frequency of photons interacting with a non- 
linear material is doubled. This phenomenon was first ob- 
served in 1961^ and has found numerous applications in 
diverse areas of science and engineering^. Physically, the 
fundamental (pump) wave propagating through a crys- 
tal with x'^^'' nonlinearity (due to the lack of inversion 
symmetry) generates a nonlinear polarization which os- 
cillates with twice the fundamental frequency radiating 
an electromagnetic field with this doubled frequency. In 
electronics, frequency doubling is a fundamental opera- 
tion for both analog and digital systems, which is however 
achieved via complex circuits made of both passive and 
active circuit elements^. 

As we demonstrate in this Letter, the possibility of 
generating frequency doubling need not be limited to op- 
tical processes in crystals, or require complex circuits. In 
fact, we show it can be realized via a completely different 
physical mechanism using the spin Hall effeclj^. Our idea 
is to use a material with inhomogeneous doping in the di- 
rection perpendicular to the main current flow. As it was 
recently shown^, a dc voltage applied to such a system 
results in a transverse voltage, similar to the Hall voltage, 
but with a different symmetry: the sign of the transverse 
voltage due to the spin Hall effect does not depend on 
the polarity of the applied field. Therefore, when an ac 
voltage is applied, the transverse voltage oscillations are 
similar in the positive and negative half-periods of the 
applied voltage, which results in the frequency doubling. 

Moreover, the transverse voltage oscillations show hys- 
teretic behavior at different frequencies. This result 
is reminiscent of the recent experimental demonstra- 
tion of memory-resistive (memristive) behavior in cer- 
tain nanoscale systems-'- and is consistent with our 
suggestion^ that some semiconductor spintronic systems 
are intrinsically memristive systems^. When a time- 
dependent voltage is applied to such systems, their re- 
sponse is delayed because the adjustment of spin po- 
larization to changing driving field requires some time 



(due to spin relaxation and diffusion processes). In other 
words, the electron spin polarization has a short-time 
memory on its previous state*. A unique feature of the 
system investigated in this work is that effects of spin 
memory manifest themselves in the voltage response, 
while in the previous study* spin memory effects were 
predicted in the current. Below, we study the frequency 
doubling and manifestation of spin memristive effects 
both analytically and numerically. This work reveals fun- 
damental aspects of the spin Hall effect that have not 
been explored yet, as well as its possible use in electronic 
circuits. 

Fig. [1] shows a possible experimental setup which can 
be used to observe frequency doubling by using the spin 
Hall effect. An alternating voltage V{t) = Vb cos {tot) of 
frequency u and amplitude Vq is applied along a sam- 
ple of semiconductor material (we call this x direction). 
The electron density in the semiconductor is inhomoge- 
neous in the direction (y) perpendicular to main current 
flow. A transverse time-dependent voltage Vt then de- 
velops between the sample boundaries y = and y — L, 
where L is the sample width. According to our previous 
calculations^ and recent experimental results^, it is ex- 
pected that the transverse voltage oscillation amplitude 




FIG. 1: (Color online) Schematics of the experimental setup. 
Inhomogeneous charge density is represented via a color gra- 
dient. 
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is a non- linear function of Vq. 

To study this system, we employ a self-consistent two- 
component drift-diffusion mode l^^'^^ which is appropri- 
ate for the physics of this problem. The inhomogeneous 
charge density profile n(y) is defined via an assigned 
positive background density profile N{y) (such as the 
one shown in Fig. [T]), which can be obtained in differ- 
ent ways including inhomogeneous doping, variation of 
sample height or gate-induced variation of electron den- 
sity. Assuming homogeneous charge and current densi- 
ties in the x direction and homogeneous z-component of 
the electric field in both x and y directions, the set of 
equations to be solved is 



Analytical solution - Before solving Eqs. ©-(O nu- 
merically, an instructive analytical result can be ob- 
tained in the specific case of exponential doping pro- 
file N{y) = Aexp(ay), with a a positive constant. At 
small values of Eq, we search for a solution in the form 
n = nQ + Sn, Ey = Ey^ + 5Ey, P ^ Pq + 5P. Setting 
Eq = 0, the leading terms in the above expansions can 
be easily obtained (see also Ref. Isj): uq = N{y) = Ae"^, 
Ey — — Po = 0. Next, using Eq. ^ and neglecting 
the term ^ ShEq, we obtain 
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where — e is the electron charge, Jt-kx) is the density of 
spin-up (spin-down) electrons, iy,|(x) is the current den- 
sity, Tsf is the spin relaxation time, crm) = en|(x)/i is 
the spin-up (spin-down) conductivity, /j, is the mobil- 
ity, D is the diffusion coefficient, e is the permittivity 
of the bulk, and 7 is the parameter describing deflection 
of spin-up (+) and spin-down (-) electrons. The cur- 
rent Ix,^{i) in a;-direction is coupled to the homogeneous 
electric field E it) = Eq cos [ujt) in the same direction 
as Ix,-\{i) = cn|(x)/i£' (i). The last term in Eq. Q is 
responsible for the spin Hall effect. 

Equation ([1]) is the continuity relation that takes into 
account spin relaxation and Eq. ([3]) is the Poisson equa- 
tion. Equation Q is the expression for the current in y 
direction which includes drift, diffusion and spin Hall ef- 
fect components. We assume here for simplicity that Z?, 
jjL, Tsf and 7 are equal for spin-up and spin-down elec- 
tronsii^ In our model, as it follows from Eq. the 
spin Hall correction to spin-up (spin-down) current (the 
last term in Eq. ^ is simply proportional to the local 
spin-up (spin-down) density. All information about the 
microscopic mechanisms for the spin Hall effect is there- 
fore lumped in the parameter 7. 

Combining equations ([T]) and ^ for different spin com- 
ponents we can get the following equations for the elec- 
tron density n — ni^ + and the spin density imbalance 
P = 7i| - ni: 
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Combining Eqs. fS]), (|4|), and integrating in y 
and neglecting the term proportional to SnSEy (this ap- 
proximation neglects small-amplitude higher harmonics 
terms), we obtain the following equation for SEy 
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Equation ([7]) already demonstrates that the driving term 
for 6Ey involves a doubled frequency. In order to further 
proceed, let us consider a sample of a finite width L, in 
which doping level variations are not dramatic. Then, in 
the first term in the right hand side of Eq. ([7]), we can 
write approximately e"^ ~ e"** , where < y* < L. This 
approximation allows us to find 

6Ey = (^Ci + \lc^ + Cl cos(2tjt - 6)^ e"*', (8) 

where 9, defined as tan 6* = C2/C3, is a phase shift, and 
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Finally, the transverse voltage is given by 
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Vt = - J SEydy 
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Equation demonstrates that there are two contri- 
butions to the transverse voltage Vt'- a shift term (pro- 
portional to Ci) and a double-frequency phase-shifted 
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FIG. 2: (Color online) Frequency doubling in transverse volt- 
age oscillations calculated at u — lO^Hz. A phase shift be- 
tween the applied field E{t) solid line) and the transverse 
field Vt (red solid line) oscillations is clearly seen. Inset: the 
transverse voltage oscillation amplitude as a function of the 
applied voltage frequency (the dots represent calculated val- 
ues of A Vt, the solid line is a fit to this points) . The plots were 
obtained using the parameter values fi — 8500cm'^/(Vs), D = 
55cmVs, e = 12.4, r,f = 10ns, 7 = 10-^ Eg = lOOV/cm 
and the background density profile A'' = 10^^ exp{2y / L)cm~'^ , 
where L = 100/xm is the sample width and < y < L. 
The analytical curve (black dashed line) was obtained at 
y* = L/1.65. 



oscillation term (proportional to \/cf"^i-C^). We have 
found that Eq. is in excellent agreement with results 
of our numerical calculations (given below) with the only 
one adjustable parameter y*. For a particular set of pa- 
rameters used below, a perfect match between analytical 
and numerical calculations was obtained at y* — L/1.65. 

Numerical solution - Equations ([3])- ([5]) can be solved 
numerically for any reasonable form of N{y). We choose 
an exponential profile for its simplicity, the possibility to 
realize it in practice, and for the purpose of comparison 
with the above analytical results. We solve these equa- 
tions iteratively, starting with the electron density n{y) 
close to N{y) and P{y) close to zero and recalculating 
Ey{y) at each time stepJ^ At each time step, the trans- 
verse voltage is calculated as a change of the electrostatic 
potential across the sample. 

Selected results of our numerical calculations are pre- 
sented in Figs. [2] and [H In particular. Fig. [2] demon- 
strates that the transverse voltage oscillations are indeed 
of a doubled-frequency character and phase-shifted with 
respect to the applied voltage. Another important fea- 
ture shown in Fig. [2] is the excellent agreement between 
our analytical and numerical results. This agreement was 
obtained by an appropriate choice of the parameter y* de- 
fined after Eq. ([7]). We observed that y* slightly depends 
on Eq and uj, and, once y* is selected, the numerical and 
analytical solutions are in a good agreement in a wide 
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FIG. 3: (Color online) The transverse voltage as a function of 
applied electric field of (a) lower and (b) higher frequencies. 



range of excitation voltage parameters. 

Multiple signatures of spin memristive behavior are 
clearly seen in Figs. [2] and [3] including, constant and 
phase shifts in Vt depicted in Fig. [2l frequency de- 
pendence of the transverse voltage oscillation amplitude 
shown in the inset of Fig. ^ and hysteresis behavior plot- 
ted in Fig. [31 All these features have a common origin: 
the adjustment of electron spin polarization to changing 
voltage takes some time. In particular, at low frequen- 
cies, we observe a small hysteresis in Fig. ^a.) because 
when the applied electric field is changed slowly (on Vt 
equilibration time scale), at each moment of time t the 
instantaneous Vt is very close to its equilibrium value ir- 
respective of the driving field E{t). At high frequencies, 
the situation is opposite: when the applied electric field 
changes very fast, the electrons "experience" an average 
(close to zero) applied electric field, resulting in a signifi- 
cantly reduced transverse voltage oscillations amplitude. 

We also note that at those moments of time when 
E{t) = 0, the transverse voltage is very close to, but not 
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FIG. 4: (Color online) Frequency doubling with the inverse 
spin Hall effect where a nonmagnetic material (green) is sand- 
wiched between two ferromagnetic leads (blue). The trans- 
verse voltage oscillation frequency depends on the relative ori- 
entation of the magnetization of the ferromagnetic contacts. 
The frequency is doubled when the direction of magnetization 
of ferromagnetic contacts is antiparallel. When the direction 
of magnetization of the ferromagnetic contacts is parallel, the 
transverse voltage frequency is equal to the applied voltage 
frequency. The oscillation amplitudes are not to scale. 



exactly, zero. This small deviation from an ideal mem- 
ristive behavior^ (predicting Vr = when E{t) — and 
related to the absence of energy storage) should be a com- 
mon feature of solid-state memristive systems operating 
at frequencies comparable to the inverse characteristic 
time of charge equilibration processes. 

We conclude by noting that another way to realize the 
frequency doubling discussed in this paper can be real- 



ized by sandwiching a nonmagnetic homogeneous mate- 
rial between two ferromagnets (see schematic in Fig. 2]). 
In this case one can employ the inverse spin Hall effect, 
in which the spin current flowing in the nonmagnetic ma- 
terial induces transverse charge current, and thus causes 
charge accumulatio n^^'^^i^''''^^i-'^^ . Since the electromotive 
force in the inverse spin Hall effect isi^ ~ Js x <?i where 
Js is the spin current along the sample and a is the 
Pauli matrix, the simultaneous change of the current di- 
rection and its spin polarization, that occurs at antipar- 
allel magnetization of ferromagnetic contacts, does not 
change the electromotive force polarity leading to trans- 
verse voltage oscillations with a doubled frequency. On 
the other hand, for parallel orientation of the spin po- 
larization of the ferromagnetic contacts, the direction of 
spin current changes within each voltage cycle but not 
the the direction of spin polarization. This thus results 
in a transverse voltage with the same frequency as the 
longitudinal one, albeit with lower amplitude (see Fig. 
111). However, since spin injection allows for much higher 
levels of electron spin polarization (of the order of several 
tens of percents) compared to the spin Hall effect (nor- 
mally, less than one percent), we expect the amplitude 
of Vt oscillations to be larger in the inverse spin Hall ef- 
fect. An external magnetic field can also be used in such 
experiments as an additional control parameter. 

Finally, the phenomena we predict can be easily veri- 
fied experimentally. They provide additional insight on 
the spin Hall effect and may find application in, e.g., ana- 
log electronics requiring frequency doubling. We thus 
hope our work will motivate experiments in this direc- 
tion. 
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